The Roles of Transmembrane Domain Helix-III during Rhodopsin Photoactivation

Background: Rhodopsin, the prototypic member of G protein-coupled receptors (GPCRs), undergoes isomerization of 11- cis-retinal to all-trans-retinal upon photoactivation. Although the basic mechanism by which rhodopsin is activated is well understood, the roles of whole transmembrane (TM) helix-III during rhodopsin photoactivation in detail are not completely clear. Principal Findings: We herein use single-cysteine mutagenesis technique to investigate conformational changes in TM helices of rhodopsin upon photoactivation. Specifically, we study changes in accessibility and reactivity of cysteine residues introduced into the TM helix-III of rhodopsin. Twenty-eight single-cysteine mutants of rhodopsin (P107C-R135C) were prepared after substitution of all natural cysteine residues (C140/C167/C185/C222/C264/C316) by alanine. The cysteine mutants were expressed in COS-1 cells and rhodopsin was purified after regeneration with 11-cis-retinal. Cysteine accessibility in these mutants was monitored by reaction with 4, 49-dithiodipyridine (4-PDS) in the dark and after illumination. Most of the mutants except for T108C, G109C, E113C, I133C, and R135C showed no reaction in the dark. Wide variation in reactivity was observed among cysteines at different positions in the sequence 108–135 after photoactivation. In particular, cysteines at position 115, 119, 121, 129, 131, 132, and 135, facing 11-cis-retinal, reacted with 4-PDS faster than neighboring amino acids. The different reaction rates of mutants with 4-PDS after photoactivation suggest that the amino acids in different positions in helix-III are exposed to aqueous environment to varying degrees. Significance: Accessibility data indicate that an aqueous/hydrophobic boundary in helix-III is near G109 and I133. The lack of reactivity in the dark and the accessibility of cysteine after photoactivation indicate an increase of water/4-PDS accessibility for certain cysteine-mutants at Helix-III during formation of Meta II. We conclude that photoactivation resulted in water-accessible at the chromophore-facing residues of Helix-III.National Institutes of Health (U.S.) (grant GM28289)National Eye Institute (Grant Grant EY11716)National Science Foundation (U.S.) (grant EIA-0225609

Entropic Tension in Crowded Membranes

Unlike their model membrane counterparts, biological membranes are richly decorated with a heterogeneous assembly of membrane proteins. These proteins are so tightly packed that their excluded area interactions can alter the free energy landscape controlling the conformational transitions suffered by such proteins. For membrane channels, this effect can alter the critical membrane tension at which they undergo a transition from a closed to an open state, and therefore influence protein function \emph{in vivo}. Despite their obvious importance, crowding phenomena in membranes are much less well studied than in the cytoplasm. Using statistical mechanics results for hard disk liquids, we show that crowding induces an entropic tension in the membrane, which influences transitions that alter the projected area and circumference of a membrane protein. As a specific case study in this effect, we consider the impact of crowding on the gating properties of bacterial mechanosensitive membrane channels, which are thought to confer osmoprotection when these cells are subjected to osmotic shock. We find that crowding can alter the gating energies by more than $2\;k_BT$ in physiological conditions, a substantial fraction of the total gating energies in some cases. Given the ubiquity of membrane crowding, the nonspecific nature of excluded volume interactions, and the fact that the function of many membrane proteins involve significant conformational changes, this specific case study highlights a general aspect in the function of membrane proteins.Comment: 20 pages (inclduing supporting information), 4 figures, to appear in PLoS Comp. Bio

Identifying Ligand Binding Conformations of the β2-Adrenergic Receptor by Using Its Agonists as Computational Probes

Recently available G-protein coupled receptor (GPCR) structures and biophysical studies suggest that the difference between the effects of various agonists and antagonists cannot be explained by single structures alone, but rather that the conformational ensembles of the proteins need to be considered. Here we use an elastic network model-guided molecular dynamics simulation protocol to generate an ensemble of conformers of a prototypical GPCR, β2-adrenergic receptor (β2AR). The resulting conformers are clustered into groups based on the conformations of the ligand binding site, and distinct conformers from each group are assessed for their binding to known agonists of β2AR. We show that the select ligands bind preferentially to different predicted conformers of β2AR, and identify a role of β2AR extracellular region as an allosteric binding site for larger drugs such as salmeterol. Thus, drugs and ligands can be used as "computational probes" to systematically identify protein conformers with likely biological significance. © 2012 Isin et al

Comparative Analysis of the Heptahelical Transmembrane Bundles of G Protein-Coupled Receptors

Background: G protein-coupled receptors represent a large family of eukaryotic membrane proteins, and are involved in almost all physiological processes in humans. Recent advances in the crystallographic study of these receptors enable a detailed comparative analysis of the commonly shared heptahelical transmembrane bundle. Systematic comparison of the bundles from a variety of receptors is indispensable for understanding not only of the structural diversification optimized for the binding of respective ligands but also of the structural conservation required for the common mechanism of activation accompanying the interaction changes among the seven helices. Methodology/Principal Findings: We have examined the bundles of 94 polypeptide chains from almost all available structures of 11 receptors, which we classified into either inactivated chain or activated chain, based on the type of bound ligand. For the inactivated chains, superposition of 200 residue bundles by secondary structure matching demonstrated that the bound ligands share a laterally limited cavity in the extracellular section of the bundle. Furthermore, a distinct feature was found for helix III of bovine rhodopsin, which might have evolved to lower its activity in the presence of 11-cis-retinal, to a level that other receptors could hardly achieve with any currently available ligands. Conclusions/Significance: Systematic analysis described here would be valuable for understanding of the rearrangement o

Role of structural dynamics at the receptor G protein interface for signal transduction

GPCRs catalyze GDP/GTP exchange in the α-subunit of heterotrimeric G proteins (Gαßγ) through displacement of the Gα C-terminal α5 helix, which directly connects the interface of the active receptor (R*) to the nucleotide binding pocket of G. Hydrogen-deuterium exchange mass spectrometry and kinetic analysis of R* catalysed G protein activation have suggested that displacement of α5 starts from an intermediate GDP bound complex (R*•GGDP). To elucidate the structural basis of receptor-catalysed displacement of α5, we modelled the structure of R*•GGDP. A flexible docking protocol yielded an intermediate R*•GGDP complex, with a similar overall arrangement as in the X-ray structure of the nucleotide free complex (R*•Gempty), however with the α5 C-terminus (GαCT) forming different polar contacts with R*. Starting molecular dynamics simulations of GαCT bound to R* in the intermediate position, we observe a screw-like motion, which restores the specific interactions of α5 with R* in R*•Gempty. The observed rotation of α5 by 60° is in line with experimental data. Reformation of hydrogen bonds, water expulsion and formation of hydrophobic interactions are driving forces of the α5 displacement. We conclude that the identified interactions between R* and G protein define a structural framework in which the α5 displacement promotes direct transmission of the signal from R* to the GDP binding pocket

Free backbone carbonyls mediate rhodopsin activation

Conserved prolines in the transmembrane helices of G-protein-coupled receptors (GPCRs) are often considered to function as hinges that divide the helix into two segments capable of independent motion. Depending on their potential to hydrogen-bond, the free C=O groups associated with these prolines can facilitate conformational flexibility, conformational switching or stabilization of the receptor structure. To address the role of conserved prolines in family A GPCRs through solid-state NMR spectroscopy, we focus on bovine rhodopsin, a GPCR in the visual receptor subfamily. The free backbone C=O groups on helices H5 and H7 stabilize the inactive rhodopsin structure through hydrogen-bonds to residues on adjacent helices. In response to light-induced isomerization of the retinal chromophore, hydrogen-bonding interactions involving these C=O groups are released, thus facilitating repacking of H5 and H7 onto the transmembrane core of the receptor. These results provide insights into the multiple structural and functional roles of prolines in membrane proteins

Site-Directed Mutations and the Polymorphic Variant Ala160Thr in the Human Thromboxane Receptor Uncover a Structural Role for Transmembrane Helix 4

The human thromboxane A2 receptor (TP), belongs to the prostanoid subfamily of Class A GPCRs and mediates vasoconstriction and promotes thrombosis on binding to thromboxane (TXA2). In Class A GPCRs, transmembrane (TM) helix 4 appears to be a hot spot for non-synonymous single nucleotide polymorphic (nsSNP) variants. Interestingly, A160T is a novel nsSNP variant with unknown structure and function. Additionally, within this helix in TP, Ala1604.53 is highly conserved as is Gly1644.57. Here we target Ala1604.53 and Gly1644.57 in the TP for detailed structure-function analysis. Amino acid replacements with smaller residues, A160S and G164A mutants, were tolerated, while bulkier beta-branched replacements, A160T and A160V showed a significant decrease in receptor expression (Bmax). The nsSNP variant A160T displayed significant agonist-independent activity (constitutive activity). Guided by molecular modeling, a series of compensatory mutations were made on TM3, in order to accommodate the bulkier replacements on TM4. The A160V/F115A double mutant showed a moderate increase in expression level compared to either A160V or F115A single mutants. Thermal activity assays showed decrease in receptor stability in the order, wild type>A160S>A160V>A160T>G164A, with G164A being the least stable. Our study reveals that Ala1604.53 and Gly1644.57 in the TP play critical structural roles in packing of TM3 and TM4 helices. Naturally occurring mutations in conjunction with site-directed replacements can serve as powerful tools in assessing the importance of regional helix-helix interactions

Helix movement is coupled to displacement of the second extracellular loop in rhodopsin activation

The second extracellular loop (EL2) of rhodopsin forms a cap over the binding site of its photoreactive 11-cis retinylidene chromophore. A crucial question has been whether EL2 forms a reversible gate that opens upon activation or acts as a rigid barrier. Distance measurements using solid-state 13C NMR spectroscopy between the retinal chromophore and the β4 strand of EL2 show that the loop is displaced from the retinal binding site upon activation, and there is a rearrangement in the hydrogen-bonding networks connecting EL2 with the extracellular ends of transmembrane helices H4, H5 and H6. NMR measurements further reveal that structural changes in EL2 are coupled to the motion of helix H5 and breaking of the ionic lock that regulates activation. These results provide a comprehensive view of how retinal isomerization triggers helix motion and activation in this prototypical G protein-coupled receptor. © 2009 Nature America, Inc. All rights reserved

Ligand-Dependent Conformations and Dynamics of the Serotonin 5-HT2A Receptor Determine Its Activation and Membrane-Driven Oligomerization Properties

From computational simulations of a serotonin 2A receptor (5-HT2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization

Understanding the mutation-induced activation of the lutropin receptor from computer simulation

The luteinizing hormone receptor (LHR) is a member of the superfamily of G protein-coupled receptors and, in humans, binds two closely related ligands, members of the heterodimer glycoprotein hormone family. This receptor is an essential component of the reproductive axis in males and females and is particularly prone to spontaneous pathogenic activating and inactivating mutations.3D-model building and computer simulations have been employed to study the mutation-induced activation mechanism of the LHR.The results of molecular simulations on different LHR models converge into the hypothesis that the arginine of E/DRY/W sequence is an important switch of the LHR activation. They also suggest that a structural modification at the interface between the cytosolic extensions of helix 3 and helix 6 is important in the mutation-induced LHR activation and/or G protein recognition. The theoretical models provide insights into the structural features of the LHR sites susceptible to spontaneous activating mutations, constituting also useful tools for "in silico" prediction of the functional behaviour of LHR mutants